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1. Introduction 

Dispersion interactions between macroscopic bodies have been the subject of 
investigation for a long time, and a variety of experiments have been performed to test 
the theoretical predictions. Whilst the theory of Casimir forces [1] between bodies in free 
space is well defined, the issue of choosing the correct stress tensor for bodies immersed in 
dielectric backgrounds is not a closed discussion [21 [3j HJ El [6] . Depending on what type of 
compatibility is required with the microscopic degrees of freedom, one is led to different 
choices. The classic example of the Minkowski vs. canonical stress tensor (formally 
the same as the stress tensor in microscopic electrodynamics), or the choice for the 
momentum density vs. Poynting vector describing momentum density propagation in a 
dispersive media show the different possibilities on offer. The canonical stress tensor is 
the choice compatible with the Lorentz force law [2] that describe the microscopic degrees 
of freedom. With this choice, the formalism of macroscopic quantum electrodynamics [7J 
can be used to evaluate the two point correlation functions. In this description one is 
able to maintain the equal time commutation relations for the physical fields and thus 
employ a sensible quantum theory of light interacting linearly and causally with matter. 
In the same fashion, many calculations are successfully carried out in the well established 
Lifshitz framework [H |9l [10] . It is also not surprising that there should be a large overlap 
with colloid physics and the corresponding dispersion forces (see for instance [llj). Given 
that Casimir force experiments can now be performed with bodies immersed in dielectric 
media [121 131 EH], it is an important objective to develop the predictions based on the 
canonical stress tensor for future experimental matching. With this testing in mind, 
we may start to think of new physical probes e.g. many body interactions [15] and 
critical Casimir effects [16J as a way of better understanding and testing the underlying 
processes. This combined effort is necessary in establishing any scheme for the choice 
of the relevant physical observables. 

The necessity of a detailed understanding of the above issues is clearly of particular 
relevance to nano-scale architectures, where a precise knowledge of these forces (stiction 
forces) between nano-mechanical objects is important. At such small length scales 
where bodies are attracting/repelling due to these interactions, one will encounter new 
frictional forces or many body effects. Atom chip miniaturisation [17] is a good example 
where one would need a quantitative description. We also see the point here that 
a detailed understanding of both different geometries and many body effects will be 
required. 

Clearly the practical need is to better understand the nature of the interactions 
between macroscopic bodies. As recently pointed out in [18] one needs to go beyond 
the simple proximity force approximation (PFA) and correctly include volume (bulk) 
contributions to the Casimir force (see also pi)] for an overview). To this end, Casimir 
energies have been calculated recently [191 ISHl [ZH [221 [23] between arbitrary compact 
bodies. Their approach was to evaluate a suitable energy functional by the use of 
path integral determinants and to make expansions, having evaluated a T-operator 
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(see also f2M 1251 [25] for similar early work). By deducing an interaction potential, 
a force could then be derived by simple differentiation (for an overview, see e.g. 
Ref. [27]). Related work on both calculational and matters of principle can be found 

in [2BII2S1 EQl ED E2]. 

In this article we evaluate the Casimir force for a collection of dielectric spheres 
in a dielectric media using the framework of macroscopic quantum electrodynamics as 
applied in Ref. [2] to the planar geometry. We use the canonical stress tensor to calculate 
the resultant force on one of the spheres in the configuration. A brief comment will be 
made on how they qualitatively differ from those that arise from the Minkowski stress 
tensor, though a detailed comparison of quantitative differences will be left to the future 
(when the appropriate experiments have the necessary sensitivity). The permittivities 
are considered to be frequency dependent and complex in general. The examples chosen 
to illustrate the inter-sphere forces will be evaluated in the retarded limit whereby we 
will use their constant static values. The calculations are carried out in a similar fashion 
to the Casimir energies cited above, by performing perturbative expansions of scattering 
two-point functions in terms of the relevant reflection/scattering coefficients. That these 
are the relevant perturbation parameters follows directly from energy conservation (i.e. 
the dielectric spheres are not amplifying bodies) and that they are effectively the energy 
dependent coupling constants of the theory. Generally speaking, one needs to find an 
approximate geometry for which the Helmholtz equation can be solved (i.e. we perform 
the separation of variables in local coordinate patches and glue them together). Then 
the eigenf unctions and eigenvalues can be used to evaluate the scattering correlation 
functions in one of two ways. One is to use a Lippmann-Schwinger evaluation of 
the Green's function based on the known background Hamiltonian where the dielectric 
spheres are perturbation potentials on top of this. The second route which is the one 
followed here is to apply the boundary conditions at all the interfaces and then deduce 
the 'out' states (the scattered modes together with their outgoing wavefunctions) from 
the 'in' states (the incident driving modes and their incoming wavefunctions). The 
driving modes, which are the coefficients of the eigenfunction expansion of the electric 
and magnetic fields (the precise form will be given later on), must be evaluated and 
indeed these are deduced from the background when no spheres are present. 

In fact, although one usually talks about the separation of the spheres, it is the 
total path length that is the meaningful distance when talking about the iV-body force. 
A further simplification is possible for a set of terms in the multiple scattering expansion 
where the number of spheres becomes large, whilst the coupling constants become small. 
For some fixed separation of the spheres we can deduce an alternative functional form 
of the force. 

The article is organised as follows. In Sec. 12] we evaluate the scattering two- 
point functions for the given setup of N distinct spheres as a perturbative series in 
the Mie scattering coefficients. This allows us to write the Casimir force as a multiple 
scattering expansion. In Sections [3] and H] we consider the two-sphere and three sphere 
configurations both at zero and finite temperatures. The evaluation performed here is 
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for the retarded limit where we take the static values of the permittivities, as a definite 
way in which to evaluate frequency integral. In Sec. [5] we calculate the iV-body force in 
the limit of a large number of spheres and at weak coupling. Finally in Sec. |6] we set 
out our conclusions. 

2. The general force expression in an N-sphere system 

It is worth remarking on some aspects of the geometry for the simple two sphere system 
before considering the general iV-sphere system. One would expect the force on one 
of the two spheres (with separation r) can be evaluated in three distinct perturbative 
regions i.e. particular distinct separations. In each region the force will have some 
scaling law:- 

(i) Region 1, where the separation of the spheres is much greater than the sphere radii. 
Here an asymptotic multipole expansion is relevant. The polarizabilities act as the 
perturbation theory coupling constants. The leading term for the force goes as 1/r 8 
and the effective geometry is three dimensional. 

(ii) Region 2, where the separation of the sphere surfaces is small compared to the radii. 
The leading term for the force per unit area goes as 1/r 4 for small parallel elements 
of area. The effective geometry becomes planar (one dimensional). This would be 
the basis for the proximity force approximation (PFA) where curvature corrections 
are then taken into account at close separations. The force then goes as 1/r 3 . 

(iii) Region 3, where the separation is comparable to the radii. Here one might think 
that there exists a third perturbative point (effective geometry) that is somewhere 
between the forms of regions one and two. The leading term for the force per 
unit length may be expected to go as 1/r 6 with the corresponding two dimensional 
effective geometry. 

The 'effective geometry' mentioned above is just the space in which the bodies can be 
moved and that parametrize the interaction. For example the parallel plate geometry 
has just the one dimension in which they can be moved. By interpolating the known 
perturbative solutions in regions 1 and 2 one can find an approximate solution that can 
be used over all distances. However, this method also requires the use of two different 
cutoff separations, where in the lower limit we would encounter an infinity and an upper 
cutoff to say what we mean by a 'large' separation. Is is also difficult to glue these 
asymptotic solutions together in the intermediate region because we need to convert the 
small separation pressure into a force. It is not obvious what area integration to use 
here and signals the trouble of using the PFA in high curvature compact geometries. 

2.1. The scattering mode decomposition 

The general strategy in calculating the Casimir force on a particular body is simply 
to calculate the scattering 2-point correlation functions for the physical fields, in our 
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case the electric and magnetic fields. Given iV-spheres we are going to calculate the 
force on sphere 1 due to the presence of the remaining (N — 1) spheres. The Helmholtz 
equation for the electric and magnetic fields is solved in N separate coordinate systems 
centred on each sphere in terms of vector wave functions about each of the separate 
spherical coordinate systems origin [33], [31]. The pictorial setup of our system is given 
in Figure [TJ 

SPHERE 2 SPHERE 3 




Figure 1. The N-Sphere system consists of N dielectric spheres of radii R[l] , • • • , R[N] 
each centred on N separate coordinate systems Si,---,Ejv, all contained in a 
background dielectric. 



Using the results in |33j, one is able to map solutions about any given sphere (any 
coordinate system Sj) into a solution about any other translated origin (a coordinate 
system E,,-). It is this fact that will enable us to perform the integral over the volume 
of the sphere. 

The Casimir force on sphere- 1 (with volume B 2 , the ball which has the two- 
dimensional sphere as its boundary) due to the effects of the N-sphere system of differing 
material properties is given by 

Fj(l\N-l)= [ d 3 xViTij(x). (1) 
Jb 2 

The stress-energy tensor is given by the standard vacuum expression (which is consistent 
with the Lorentz force law [2]) 

Tijix) = E i (x)E i (x) + B i (x)B i (x) - ~<%(|E(x)| 2 + |B(x)| 2 ), (2) 
and it is understood that we are taking the limit for the initial and final points 

lim E(a; 1 )E(x 2 ) = E(ar 2 )E(ar 2 ), (3) 

lim B(x l )B(x 2 ) = B(x 2 )B(x 2 ). (4) 

Il->12 

In performing this limit it is first necessary to make subtractions of the singular part 
arising from the free Green's tensor. This is somewhat like a vacuum bubble one 
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encounters in quantum field theory where we are throwing away an unobservable zero 
point energy. We then need to evaluate the scattering correlation functions (whilst 
dropping the direct modes of propagation), viz. 

POO 

limE i (x)E J -(y)= / du>du>\E^(x;uyE i ?(y;u/)), (5) 

V-** Jo 

and similarly for magnetic fields. 

A crucial second step is to realise that the volume integral of the force requires a 
careful specification of the coordinate origin in order to be able to evaluate the derivative. 
In the multiple scattering approach used here the coordinate origin used will be the one 
centred on the sphere from which the last scattering event took place, then translated 
to the origin of the sphere where we are calculating the force. In terms of the classical 
scattering Green's tensor Gy centred on the sphere where we are calculating the force, 
schematically we have 

/ d 3 sV i G ii (x,x) « R I dfi 2 i? 2 V x G rr (R,R), (6) 
Jb 2 Js 2 

where X is the vector connecting the centre of the two spheres, originating from the last 
scattering sphere centre and dVl 2 is the element of solid angle. It should be pointed out 
here that this approximation, in order to be consistent with the case of a sphere and a 
dipole, requires the introduction of a normalisation factor (for us it will be l/Air). This 
is because the derivative acting on the stress tensor is evaluated at the centre of the 
sphere, whilst the wavefunctions are evaluated on the spheres surface. In the case of the 
sphere and dipole, we do not have the integral to do over the S 2 . Due to the potentially 
singular nature when shrinking the sphere down to a point, this clearly must be done 
carefully. 

To construct the scattering two point function we write the fields in a mode 
decomposition |33j of spherical vector wave functions 

|E m ) = Yl Pi, m ML,m,j L ) + £ €m A L) \L,mJ L ), (7) 

i,L,m i,L,m ^ 

|E out )= J2<,mML,m,ht)+ E 6 UQ VAL ) l^ht), (8) 

i,L,m i,L,m ^ ' 

|E""> = A >m ML,m,j L ) + 4, m AL) \L,m,j L ), (9) 

i,L,m i,L,m 

|B) = -VA|E). (10) 

(Note that the mode operator coefficients are functions of frequency which for now has 
been suppressed). The 'in' and 'internal' states are regular at the i-sphere origin, whilst 
the 'out' states are outgoing modes falling off at infinity. They are eigenfunction modes 
with respect to the i-sphere of the radial eigenfunctions, represented by spherical Bessel 
and Hankel functions respectively. The internal coefficients (c l Lm , d l L m ) of a given sphere 
can be eliminated by application of the continuity equations across each sphere. One 
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is then able to determine the scattering coefficients (a l L m , b l L m ) in terms of the driving 
modes {PLmilLm) by making use of a translation addition theorem for vector wave 
functions [33]. These are given by 



L\H,Ll,ml,j L1 )= Yl A Li% ; L2,m2 L \ i2 , L2 , m2 , h L2) 
i2,L2,m2 

+ ^ 1 i';L;L2, m2 ^VALK2,L2,m2 ) / i + 2 ), (11) 

and 

^VAL|il,Ll,ml,j L i>= ^ 1 ili;L 2im2 ^VAL|z2,L2,m2,/i+ 2 ) 

i2,L2,m2 ' 

+ Bl\% ]L2 , m2 L\t2,L2,m2,h + L2 ). (12) 

Here, A and B are representations of the translation group in the angular momentum 
basis. They can be found in [33] (see also [Mj), together with the other necessary 
quantities used here, but in particular they are functions only of the centre-to-centre 
sphere separation r[i,j] and the background k. The matrices A map TE-to-TE and TM- 
to-TM, whilst the matrices B map TE-to-TM and vice versa. Using the Mie scattering 
coefficients (ql, Pl) the scattering modes are given by (where we suppress indices) 

(a,b)(l + T) = (a-p,/3-q), (13) 

where the T-matrix has the block diagonal form 

a ■ A a ■ B 
p-A P-B 

The matrix equation given in (1131) can be inverted and then perturbatively expanded 
as a power series in the Mie scattering coefficients (which play the role of the small 
coupling constants) such that the scattered modes are found to be 

oo n 

M>=^n(- T )">-^-<?>- ( is ) 

n=0 m=0 

The T-matrix is then the set of reflection coefficients from each sphere appropriately 
translated and the integer n gives the number of scattering events for the whole set of 
spheres. This is the multiple scattering approach (see for example [33]) which allows for 
all possible scattering between the spheres. 

In order to calculate the iV-body Casimir force it is necessary to consider all possible 
ways of connecting all of the spheres. We can form a primitive set of diagrams that can 
be built on to generate all other diagrams. A suitable set of definitions are:- 

• Disconnected diagrams are where some subset of the spheres do not take part in 
the multiple iV-body scattering process. 

• Simply-connected diagrams are where all the spheres are connected with no more 
than two lines attached to any sphere. This forms a continuous path through all 
the spheres. 



(14) 
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• Reducible diagrams are where diagrams can be written as simply connected diagram 
convoluted with a lower body diagram. This affords a description of multiple 
internal scattering. 

The simply-connected diagrams can be seen to be irreducible iV-body diagrams. Unless 
otherwise stated, we will always be calculating the simply-connected diagrams and 
corresponds to one particular perturbation series (as opposed to a perturbation series 
based on the sphere radii). 



2.2. Dielectric backgrounds 

We now specialise to the case of dielectric backgrounds and consider an approximation 
where (3l = so that the b modes vanish. This approximation is based on the eventual 
small argument expansion of the Mie scattering coefficients (and thereby a simple 
multipole expansion) where the scattering coefficient gives the leading contribution 
in the frequency integral. What is interesting in this case is that there are only 
contributions from the cross terms involving the a and p modes, as the q and a cross 
terms are zero due to the antisymmetry properties of the translation products A n ■ B. 
From Equation (Tl5|) it is possible to evaluate the 2-point functions in terms of the driving 
modes on the surface of sphere-1. The relevant components (suppressing SO(3) indices) 
contributing to the rr component of the stress tensor are:- 

Js 2 

J ^^(x^B^^u/))^ = 3^ h(k\Ri\)ht(k\Rx\ 

[ (E n >™\x;uy .E i2 ' in {x;u'))\ x=m = $L(L + l)^^^^^!) 
Js 2 

■(a^y-^f-pfiu')), 

[ (B^^coY ■ B^ n (x-u'))\ x=R[1] = %- 2 (L(L + tyhWRtDh+ikW) 
Js 2 u 

Here we have used the two relations that the eigenfunctions in the bulk background 
satisfy:- 

L\L,m,j L (k\x\)) =j L (k\x\)L\L,m), 

1^ T \ lr / 1 1 i \ \ L(L + l)j L (k\x\) lT , 

-VAL L,m,jx(feN)> = — ,|| \L,m). 

k J r k\x\ 

In addition we are anticipating a ^-function in the frequency so that we have been a 
little cavalier with the frequencies above. 
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2.3. Evaluating the background vacuum 

What remains to be calculated are the expectation values of the p-p vacuum modes, 
(0|(PLimi)(Pi2m2)*|0) w ^ n respect to some sphere i (in our case this will be sphere-1). 
This can be evaluated for a dielectric background €b(u>) that is independent of position 
and when there are no spheres present. We can then perform the integral over all space. 
Note this is different from the case of two static dipoles interacting as it is only the 
location of the dipole source where the noise current is non-zero as opposed to here 
where the entire background is filled with noise which the spheres polarizability couples 
to. 

The composition rule for two free (or bulk) Green's functions in the absence of any 
perturbing spheres integrated over all space is given by (a more general result using 
functional differentiation and a proper tensor treatment is given in Appendix A ) 

d 3 zG^ e (x,z\k)(G^ e (y,z\k)y = ±-d k Gf ree (x,y\k), (16) 

with k = \J £b{u)uj / c. Using the explicit form of the free propagator for incoming modes 



i e ik\x—y\ 

'ij^{x,y\k) = [dij + ■ ' 

which becomes 



G^(x,y\k) = ^ + ^ j )-^— {} M7) 



- id k G{; ee (x,y\k) = + iViV,^. (18) 
This can then be decomposed as a spherical wave expansion (see e.g. [35] ) 

- id k G{; ee {x,y\k) = M*ij + ^2 V * V i) 

■^2jLi(k\x\ )j L2 (k\ y \)Y L1>ml Yl 2m2 . 

We can use this result to evaluate the background media vacuum that has the quantum 
noise fluctuations throughout. The in field modes are given by (at some fixed radial 
distance and in Sj coordinate system) 

PkmiM = —TTT^ I I d3zLY Lmi ■ G free (x, z\k) ■ 3(z; uj). (19) 

]L1{\X\) Js 2 J 

The quantum noise current correlation functions at zero temperature are (see [2]) 

(3l(x,u)J J (y,u')) = 0, (20) 

and 

(J^c^jW)) = - ( t 4y/$S(e B (u>))y/S(e B (u'))) 5\x - y)6(u - u'). (21) 



7T \ C 2 



After some algebra and taking the limit x — > y we obtain for the modes incident on the 
i-th sphere 

hu? 

(0\(pkml)(pkm2^\0) = ^(ik)S(u - U/)1* ® 1* , (22) 
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where 1 is the (2L + 1) vector populated entirely by l's 1 = (1, •••,].), which we 
shall from now on denote as |1). We see that there are no correlations between the 
different mode numbers as one would expect from random noise. By rotating to the 
imaginary frequency axis it is found to be a sensible positive definite object representing 
the different occupation numbers of the noise excitations. Further, the size of these 
vectors give us the level at which we are truncating our series representation in the 
angular momentum basis. 



2.4- Perturbative force for dielectrics 

By making the appropriate substitutions from the previous sections, one is able to 
calculate explicitly the force on sphere- 1 due to simply connected scattering processes 
between the iV-spheres (multiple reflections can be included in an obvious fashion). This 



is given in Appendix B Here we state simply the result: 

F [l|iV - 1] = _(_l)^ Y Vr[i y I" dX ^ : COt ( , ^ \ 

4vrf- [ ' ] Jo v^lM] \k B TD[l,l]J 



X m ■ ( XR[1] \ ( XR[1] \ ( XR[1] \ 
V[1,1] J \V[1,1]J \V[1,1]J \V[1,1]J 



|1).(23) 



Here, X>[1, 1] is the loop distance of the simply-connected path from sphere- 1 back to 
sphere-1 consisting of the sum of the N separations, 

V[l, 1] = (\r[l,i)\ + \r[i, j)\ + ■ ■ • \r[j, k)\ + \v[k, 1]|)_ 

and A l+1,1 are the set of polynomials in (X,r[i + l,z]) that remain after extracting 
the exponentials from the translation coefficients. Due to the additive nature of these 
exponentials one finds it is the total path length that occurs in the exponential. The 
function W is a measure due to the two separate contributions in the stress tensor 



arising from the individual field components and isotropic part (see Appendix B ). 

This expression can be used as a perturbative expansion. The arguments of the Mie 
scattering coefficients are always small since the total path length will always be larger 
than the sphere radii. It can then be evaluated in different perturbative setups when 
T = 0, or by residues when T ^ 0. It is also worth commenting how the expression 
for the force derived here will differ or coincide with that found from the Minkowski 
stress tensor or from energy functional methods. At least in some perturbation scheme 
there will be two main differences. One will be the scale of the force due to the differing 
factor of the background permittivity tensor. The second will be a different set curvature 
corrected terms that result from the extra isotropic contribution to the stress tensor. 
In the limit where the background permittivity is set to the vacuum value, all the 
expressions for the force will coincide in an exactly analogous way to that of the parallel 
plate geometry found in [2]. 
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3. The Force Between Two Spheres 

We now analyse the specific situation of the force on one sphere due to multiple 
scattering with second in the retarded limit. In this manner we have an illustrative 
scheme in which to evaluate the force and thereby derive simple expansions. This was 
performed recently using functional determinants in [T9], [20] where in particular the 
scattering coefficients are expanded to higher powers in the frequency (and thereby 
sphere radii). For simplicity we consider identical spheres that are aligned along the 
z-axis resulting in the diagonalisation of the m-indices and take the static values of 
the different permittivities (since we are in the retarded limit). The Mie scattering 
coefficients are also expanded for small arguments only to first order so that we will find 
static polarizabilities of the simple multi-poles entering expressions. This is done to show 
the applicability of the method. All evaluations (integration and matrix operations) have 
been performed in Maple and numerical data has been taken from [36] . 



3.1. The retarded limit in the vacuum atT = 

For the case of two spheres in the vacuum, one is able to perform an expansion in 
either the curvature of the spheres or the difference in permittivities (or the product 
combination thereof occurring in the scattering coefficients). In addition the number 
of scattering events also acts as perturbation parameter. The first perturbation series 
we consider is the difference in permittivity to make contact with known results. These 
can be assembled into the standard multipole moments plus higher curvature corrected 
contributions that are valid results at bilinear order in the difference in the permittivities. 
This second set of terms arise from the isotropic part of the stress tensor. Indeed, an 
expression is obtained directly for the two body force 

F [1121 = - V - (v a 1 a 2 + w m,n«R{i? \ m 



Here the coefficients (v m>n , w min ) are a set of numerical coefficients found after evaluating 
the frequency integral, whilst the multipole moments a l L are the multipole moments of 
sphere-i given by 

ai = e ^ e f ~ 1 ff[jfW, (25) 
e*/e B + (L+l)/L U V ; 

It should be said that the first term in the series coincides with the well known vacuum 
result with v\ t \ = 23. This has been evaluated for polystyrene balls in the vacuum with 
values given in Table [TJ The force has been plotted against separation in Figure [2j 



3.2. Finite temperature T > 

Here we consider the polystyrene balls suspended in a liquid silicone solution at room 
temperature. The values are given in Table [2l followed by their plot in Figure O The 
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Table 1. Material properties of the two sphere system. The dielectric spheres are 
polystyrene balls in a background of the vacuum at T = OK. 



Properties 


Sphere 1 


Sphere 2 


Permittivity (Sphere) 


2.6 


2.6 


Permittivity (Background) 


1.0 


1.0 


Radii (m) 


1.0 x 10~ 6 


1.0 x HT 6 



L(truncation) = 3 







- 10°- 






- 10 U 






- 10 2 -^ 
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- 10 3j 
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- 10 +J 
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- 






- lO 7 -^ 






- 10 SJ 




Inter-Sphere separation 
(in dimensionless units x = r[1 ,2]/R) 



Figure 2. A plot of the inter-sphere retarded force between two identical dielectric 
spheres (with radii R = R[l]) in the empty vacuum at T = 0°K. The vertical axis 
is the dimensionless ratio Force / (he/ (AttR 2 ) x 10 -12 ), whilst the horizontal axis is 
the dimensionless ratio x — r[l,2]/R. Here, the relative dielectric permittivity of the 
polystyrene spheres is ei = ti — 2.6 . The angular momentum series representation 
has been truncated at L trunca te = 3. 



poles are located on the imaginary frequency axis at HcXi/(kBTV[l, 1]) = In for integer 
I. The attractive force can hence be evaluated to be of the power series form 



oo oo 



x a 



1=1 m,n=l 
2m+2n 



t=0 

x ( vW + vW 



e -2r[l,2]^K B TlTr/(hc) 

i m i\ 2m+2n—t 
Kb± 7TL \ 



1 



he 

k B TR[l] \ 

he ) 



9t 



r[l,2] 



3+t 



(26) 



Once again the integration of the frequency gives another set of multipole type 
coefficients {v ; mn ,w' mn ) together with another set g t that arise due to the different 
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Table 2. Material properties of the two sphere system. The dielectric spheres are 
polystyrene balls in a background of liquid silicone at T — 290K. 



Properties 


Sphere 1 


Sphere 2 


Permittivity (Sphere) 


2.6 


2.6 


Permittivity (Background) 


2.2 


2.2 


Radii (m) 


1.0 x 1(T 6 


1.0 x KT 6 



L(truncation) = 3 




-] — i — | — . — | — i — | — i — | — i — | — i — | — i — | — 
4 6 8 10 12 14 16 18 



Inter-Sphere separation 
(in dimensionless units x = r[1,2]/R) 

Figure 3. A plot of the inter-sphere force between two identical dielectric spheres 
(with radii R = R[l]) in a silicone fluid background (e B = 2.2) at T = 293°^. The 
vertical axis is the dimensionless ratio Force/ (he/ (AttR 2 ) x 10~ 12 ), whilst the horizontal 
axis is the dimensionless ratio x = r[l, 2]/R. Here, the relative dielectric permittivity of 
the polystyrene spheres is e\ = — 2.6 . The angular momentum series representation 
has been truncated at L trunC ate = 3 and the Matsubara frequencies truncated at I = 2. 



4. The Force Between Three Spheres 

We now analyse the specific situation of the force on one sphere due to multiple scattering 
with two other spheres, again in the retarded limit for illustrative purposes. See 
also [371 138] for details of three atom potentials. It is worth commenting that in order for 
the simply connected diagrams to hold as valid, we require a clear line of sight between 
the three spheres. This means in particular that the multipole expansion will not be 
accurate when the three spheres are collinear (indeed to form a connected diagram now 
requires multiple internal reflections effectively convoluting two diagrams that are both 
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two-body). However, as with the approximations made for the translation coefficients 
(centre-to-centre neglecting the material properties of the spheres), we neglect this 
effect and assume this collinearity effect not to be present. As with the two sphere 
case, static permittivities are used (corresponding to the retarded limit) and the Mie 
scattering coefficients are expanded for small arguments so that expressions contain 
static polarizabilities. 




Figure 4. A plot of the inter-sphere retarded potential between three identical 
dielectric spheres (with radii R = R[l]) in the empty vacuum at T = 0°K. The 
vertical axis is the dimensionless ratio Potential / (Hc/(4ttR) x 1CP 11 ), the horizontal 
axis labelled by x is the dimensionless ratio x = r[l,3]/R, and the remaining axis is 
the angular variable 9[1, 3]. Here, the relative dielectric permittivity of the polystyrene 
spheres is t\ = 62 = £3 = 2.6. Spheres 1 and 2 are held fixed along the z axis 10R apart. 
The angular momentum series representation has been truncated at L trU ncate — 1- 

Figure H] shows a plot of the inter-sphere retarded potential from which the force 
can be found from simple gradients given in Equation (|23|) . This is given as a function 
of the separation and angle for the third sphere for a fixed separation vector between 
spheres 1 and 2 in the retarded limit in the vacuum at T = 0. For the finite temperature 
case T > 0, Figure |5] is a plot of the three sphere potential for a fixed separation of 
sphere 1 and 2 along the z-axis. 

5. The Large iV limit of the Inter-Sphere Force 

Consider the case of N identical spheres arranged in some random fashion as iV — > 00, 
a % = a — > 0, whilst A := Nets is held fixed (a = asUJ 3 R 3 /c 3 ). We put one constraint on 
the separation of all the spheres, that every sphere has two nearest neighbour spheres 
with a separation s between them. In this case there is a single irreducibly connected 
diagram that has a minimum path length T>[1, l] m . There are also approximately iV 
irreducibly connected diagrams that have a slightly increased path length V[l, l] m +i 
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e[i,3] 

Figure 5. A plot of the inter-sphere retarded potential between three identical 
dielectric spheres (with radii R = R[i\) in silicone fluid with eb = 2.2 at T = 293°K. 
Spheres 1 and 2 are held fixed along the z axis 10R apart. The vertical axis is the 
dimensionless ratio Potential / (Hc/(4ttR) x 10~ 20 ), the horizontal axis labelled by 
x is the dimensionless ratio x = r[l,3]/R, and the remaining axis is the angular 
variable #[1,3]. Here, the relative dielectric permittivity of the polystyrene spheres is 
e\ = 62 = £3 = 2.6 . The angular momentum series representation has been truncated 

at Ltruncate 1- 



that satisfy T>[1, l] m < V[l, l] m+ i < T>[1, l] m + 2s. We make the further approximation 
that each term of the second set contribute the same amount as the minimum path 
length diagram to the iV-body scattering. The expression for the associated potential 
reduces to 

m]« ± (-l)^ ra x* 

*f « e Mi^>^ s) i N - (27) 

Taking the T — > limit and then iV — > oo we can extract the dominant L = 1 term 

N 



(28) 



which reduces to (using Sterling's approximation) 

p -N r>3N 

V[N] « ±(-lfHc e —\ N ^. (29) 

The use of this formula could be in adding an additional sphere to the system and 
measuring the resulting oscillation of the force compared to the absolute force before 
the sphere is added. Alternatively it could serve to quantify the errors or to rearrange 
the perturbation series. 
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6. Conclusions 

We have calculated Casimir forces between spheres using a multiple scattering approach. 
This has been done at both zero and finite temperature. In the case of the simply- 
connected three-body force, we have pointed out the need for non-collinearity, but that 
it is perturbatively not necessary. The total closed path length being always greater than 
the respective length scales of the associated radii helps in quantifying a perturbative 
evaluation, for example the frequencies (found from the radii and total path length) 
below which the Mie scattering coefficients can be expanded. This is simply due to the 
fact that the propagator connecting all scattering centres in this simple setup contains 
an exponential factor of the total path length. Additionally it addresses the technical 
problem of going beyond the proximity force approximation (PFA) [18J since, although 
continuity equations are applied at each material surface, there exists a solution of the 
field equations in each material medium. In general there will be both surface modes 
and bulk excitations in the bodies and calculations of Casimir forces should contain 
both [39]. The multipole expansions found when the Mie scattering coefficients are 
expanded are also not singular when the spheres touch, since the sum of the two radii 
act as a physical lower cutoff (non-overlapping potentials). 

For large N and weakly scattering spheres, we have deduced a scaling type formula 
for the simply-connected force terms based on a fixed separation set of properties for 
the configuration. This could be useful for forces where we want to understand the 
behaviour as a function of N rather than the detailed configuration (separations and 
orientations). It can also serve to quantify the errors in a multiple scattering expansion 
when truncating the series at this order. In a subsequent paper one of us will develop 
the idea of the total path length playing a key role in Casimir interactions. 

Theoretically and experimentally there are many open issues that need to be 
addressed in conjunction and naturally follow on. It is clear that substantial progress 
has been made both theoretically and experimentally in refining the accuracy of Casimir 
force calculations and measurements [40J. However, it should be pointed out that the 
theoretical calculations are for idealised geometries that do not necessarily reflect the 
experimental reality. One is faced with a potential disparity between the geometry 
encountered in the experiment, and those computable theoretically. We may for instance 
start off with a body that is very spherical, but that when we process it (heating, gluing, 
or coating) the sphericity of it will be reduced. For instance in [13] the initial borosilicate 
spheres are extremely smooth (roughness < 0.1%) but that the process of heating can 
induce surface roughness. And even if we calculate scattering coefficients for a given 
surface roughness profile, this is not necessarily what will be found experimentally. The 
problem requires both theoretical and experimental studies that are developed together. 
It is clear that we must perform a matching of experimental data (classical scattering 
data from the body) with calculations that allow for these unknown 'to be measured' 
parameters (for example non-spherical surfaces). Only with this interplay can one be 
confident that the real bodies are being quantitatively described and to what extent 
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perfect geometries are a good approximation. It could be labelled under the 'realistic 
geometries' motif. For bodies in dispersive media, more needs to be understood on 
different material choices together with improvements experimentally, to be able to 
resolve the detailed difference between different proposals for physical observables. 
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Appendix A. Composition of Green's tensors 

The bulk Green's tensor is defined by 

(V y + kl(x,u))G{j ee (x,y\k) = 5 tJ 5 3 (x - y), (A.l) 
which can be rewritten as 

G(r(x,y\k) = (x\— -^-rrly). (A.2) 

V ij -+- K {j 

This can be functionally differentiated with respect to k^ 

I ^4^ fr(x ' m = ~ {xl (v y + m k]m '(v™Wj |y> - (A - 3) 

A complete set of states can now be inserted and together with taking the imaginary 
part we find 

/ »(*L)^-G y (x,y|fc) = - J d 3 zGU^^m(k 2 mn )(Gn,(y,z\k))l(AA) 

At this point further simplification is possible when a known Green's tensor is available 
together with the corresponding dispersion relation. 

Appendix B. Derivation of the multi-sphere force 

Substituting the scattering coefficients and the driving mode vacuum expectation values 
into the field correlation functions one finds the perturbative force on sphere 1 due to 
multiple scattering is, 



F[1|JV- 1] = - 




duk 



x (l^iuRilj^A^iv^^.a^Rli})... 

i=2 
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JV 



N 



(B.l) 



(B.2) 



...^.l-(r7.. / ;).n^V 1 ., l .K l (r,.r) 

x j(kR[l})h + (kR[l})W(coR[l])\l}, 
where (here the background permittivity is e#) 

= L(L + 1) (l(L + 1) - ^(1 + e B )u 2 R[l] 2 ^j . 

Here we have omitted the SO (3) indices as in the earlier text, but for the measure W 
above, have stated how it depends on the labels. This encapsulates all contributions 
of the force as an order N polynomial in the scattering coefficients. We use the reality 
properties of the polarizabilities and translation coefficients after Wick rotation together 
with the necessary thermal factor for the finite temperature system to rewrite this 
expression. The polynomials in the translation coefficients can be seen to be a path 
ordered expression linking all the scattering centres in a closed loop starting and finishing 
on sphere 1. We define this polynomial to be VO[a, A]^ and that it is simply connected 
(i.e. a scattering coefficient for any given sphere only appears once). With the gradient 
defined always to act on the last translation matrix, we obtain for the force, 



F[1|7V- 1] = 



{-1) N+1 H 



N 



4:71 



r[i,l] 



dQ cot 



Ml 
k;7T 



x (II 



[VO[a,A]^} 

x (kR[l]) ■ j(R[l]) ■ n(R[l])W(QR[l])\l). (B.3) 

We now make a change of integration variables to X = y/eBQT>[l, l]/c, which is 
dimensionless, together with the introduction of the Jacobian factor J and the extraction 
of the exponentials from the translation coefficients, so that 



N 



F[1|JV- 1] = - 



' =2 



-X 



X 



x (1 



X 



dxr 



■ cot 



HcX 



V[l,l] \k B TV[l,l], 
[VO[a,A]? tl ] 
[ XR[1] \ , XR[1]\ JXR[1]\ 

\v[iA\)- J \v[iA\)- \v[i,i]J 



w 



with 



J :-- 



dX 

In' 



(B.4) 



(B.5) 



For the calculations in this paper we will always take this Jacobian to be the simple 
case of a constant dielectric. 
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